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ABSTRACT 

Aims. To understand the environment and extended structure of the host galactic gas whose molecular absorption line 
chemistry, we previously observed along the microscopic line of sight to the blazar/radiocontinuum source NRAO150 
{aka B0355+508). 

Methods. We used the IRAM 30m Telescope and Plateau de Bure Interferometer to make two series of images of the 
host gas: i) 22.5" resolution single-dish maps of 12 CO J=l-0 and 2-1 emission over a 220"by 220" field; ii) a hybrid 
(interferometer+singledish) aperture synthesis mosaic of 12 CO J=l-0 emission at 5.8" resolution over a 90"-diamctcr 
region.** 

Results. At 22.5" resolution, the CO J=l-0 emission toward NRAO150 is 30-100% brighter at some velocities than seen 
previously with 1' resolution, and there are some modest systematic velocity gradients over the 220" field. Of the five 
CO components seen in the absorption spectra, the weakest ones are absent in emission toward NRAO150 but appear 
more strongly at the edges of the region mapped in emission. The overall spatial variations in the strongly emitting 
gas have Poisson statistics with rms fluctuations about equal to the mean emission level in the line wings and much 
of the line cores. The J=2-l/J=l-0 line ratios calculated pixel-by-pixel cluster around 0.7. At 6" resolution, disparity 
between the absorption and emission profiles of the stronger components has been largely ameliorated. The 12 CO J=l-0 
emission exhibits i) remarkably bright peaks, T m b = 12 — 13 K, even as 4" from NRAO150; ii) smaller relative levels 
of spatial fluctuation in the line cores, but a very broad range of possible intensities at every velocity; and iii) striking 
kinematics whereby the monotonic velocity shifts and supersonically broadened lines in 22.5" spectra are decomposed 
into much stronger velocity gradients and abrupt velocity reversals of intense but narrow, probably subsonic, line cores. 
Conclusions. CO components that are observed in absorption at a moderate optical depth (0.5) and are undetected 
in emission at 1' resolution toward NRAO 150 remain undetected at 6" resolution. This implies that they are not a 
previously-hidden large-scale molecular component revealed in absorption, but they do highlight the robustness of the 
chemistry into regions where the density and column density are too low to produce much rotational excitation, even 
in CO. Bright CO lines around NRAO150 most probably reflect the variation of a chemical process, i.e. the C + -CO 
conversion. However, the ultimate cause of the variations of this chemical process in such a limited field of view remains 
uncertain. 
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1. Introduction 

Many detailed and sensitive studies of interstellar gas em- 
ploy absorption spectra taken against the microscopic disks 
of stars or compact radio continuum sources. We have re- 
cently used such spectra at radiofrequencies to elucidate 
an unexpectedly widespr ead molecular chem istry in diffuse 
clouds; for instance, see iLiszt et ail (|2006l) and the many 
references therein). 

Lamentably, observing gas in absorption along such sin- 
gle, isolated lines of sight does not generally allow the iden- 

Send offprint requests to: pety@iram.fr 

* Based on observations obtained with the IRAM Plateau de 
Bure interferometer and 30 m telescope. IRAM is supported by 
INSU/CNRS (France), MPG (Germany), and IGN (Spain). 

** The three spectra cubes (RA, Dec, Velocity) are available in 
electronic form (FITS format) at the CDS. 



tification of the intervening features with recognizable phys- 
ical bodies in space, and their possible association with 
other nearby entities (for instance dark clouds) remains 
hidden. These limitations have obvious consequences for 
our understanding of the intervening medium and the ab- 
sorption profiles themselves. But in some rare cases it is 
becoming possible to observe the same feature in emission 
and absorption, even in the same chemical species, and so 
perhaps to characterize the nature of the intervening ma- 
terial by imaging it across the line of sight away from the 
background object. 

This paper begins a short series of reports of studies 
undertaken to elucidate the nature of the gas occulting the 
compact extragalactic mm-wave background sources em- 
ployed in our work. It describes one case, using the J=l-0 
and 2-1 mm-wave rotational transitions of carbon monoxide 
to image the particularly complex kinematics of relatively 
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distant absorbing mate rial along the l ow-latitude line of 
sight toward NRAO150 (|Cox et al J 19881) at the highest use- 
fully attainable resolution. It is organized as follows. Section 
2 describes the sightline toward the blazar NRAO150 as 
it is understood from radiofrequency observations of CO 
and HCO + in emission and in approximately a dozen other 
atomic and molecular species in absorption. Section 3 de- 
scribes the techniques used to image and/or synthesize the 
pattern of foreground CO emission using singledish and in- 
terferometeric synthesis data. Section 4 describes the re- 
sults of the new observations and Section 5 is dedicated to 
models and Section 6 is a summary. 

2. The line of sight toward the blazar NRAO150 

Unlike the vast majority of strong radio continuum 
sources, no visual identifi cation of the blazar NRAO150 
(|Paulinv-Toth et"aT] I1966D has yet been made. At a = 
3h59m29.73s, S = 50°57'50.2" (J2000), NRAO150 is seen at 
low latitude in the outer Galaxy (I, b =150.3772°, -1.6037°), 
so the path through the galactic disk is long and there is 
a heavy accumulation of distributed diffuse/H I gas. Quite 
aside from any dense material which might intervene, the 
integrated H I emission of th e nearest spectrum in the 
Leiden-Dwingeloo H I survey (|Hartmann fc Burton! Il997l ) 
gives an optically-thin lower limit N(H I) > 7.5 x 10 21 cm~ 2 
or A v > 4 mag. There is an approximate 5 kpc path length 
from the Sun out to a galactocentric radius of 13 kpc where 
the neutral gas disk might end (assuming Rq = 8.5kpc). 
For a typical mean gas density of 1 cm~ 3 , this corresponds 
to N(H) = 10 22 cm" 2 or A v w 6 mag. 

The sightline also harbors denser H2-bearing clouds 
as shown in Fig. [T] While the H I absorp tion spectrum 
is to o heavily saturated to be revealing (|Liszt fe Lucasl 
|1996|) . the CO and HCO + absorption spectra reveal 5- 
6 absorption features at LSR velocities of —17.5, —14, 
— 10.5, —8.5, —4.5/ — 3.5kms _1 . All those absorption fea- 
tures have comparable N(HCO + ) w 1 x 10 12 cm" 2 im- 
plyin g N(H) ft 2N(H 2 ) ft 2 x 10 21 cm" 2 (jLucas fc Lisztl 
1996; iLiszt fc Lucasl 119960 or A v w 1 mag. Each cloud 
also has N(OH) comparable to that seen toward ( 
Oph (ibid), t he prototypical diffuse cloud sightline stud- 
ied optically (|Van Dishoeck fc Blackl |1986[ ) where the vi- 
sual extinct ion is 1 mag. As in optically-studied dif- 
fuse clouds (jSonnentrucker et alj l2007t iBurgh et all l2007t 
iSheffer et al.ll2007l). w e derived NCCO) < 10 15 - 10 16 cm" 2 
(jLucas fc Lisztl I1998D . This implies that only a modest 
fraction of the expected free gas- phase carbon N (C) ft 
1.6 x 10" 4 N(H) ft 3 x 10 17 cm" 2 (jSofia et al.ll200l is in- 
corporated in CO, as in any cloud with Ay^ 1 mag. 

We conclude that the molecular features seen toward 
NRAO150 are hosted in diffuse or marginally translucent 
gas, consistent with the weakness of emission spectra in 
species other than CO and H I. OH emission seen at 18' 
resolution is quite weak (0.01 K) and very broad, extend- 
ing over — 35 < v < Okms" 1 and individ ual cloud features 
cannot be discerned (jLiszt fc Lucasl 1996). The only species 
seen in mm-wave emission beside CO is HCO + which ap- 
peared very weakly (0.01-0.02 K) at v = — 9kms -1 when 
avera ged over a several arcminute region (jLucas fc Lisztl 
[199(1 . 

The absorption features display a rich and varied chem- 
istry, a small sample of which is illustrated in the spec- 
tra of the bottom panel of Fig. [TJ They harbor a promi- 
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Fig. 1. Emission and absorption spectra toward 
B0355+508 aka NRAO150. Top panel: 1' HPBW 
J=l-0 12 CO and 13 CO profiles from the former NRAO 
12m telescope on a Tr* brightness scale and the 22.5"- 
resolution IRAM 30m 12 CO profile from the present work 
on a main-beam brightness temperature scale. Bottom: 
IRAM PdBI absorption profiles from various representative 
species. The profiles have been shifted vertically and scaled 
in some instances. 



nent complement of typically ubiquitous species such as 
CO, C 2 H and 3h2 (jLucas fc Lisztl l200l as well as OH. 
All of the clouds also have detectable absorption in HNC, 
although weakly at — 14kms _1 and —4.5/ — 3.5kms _1 . 
The yet-stronger HCN spectrum is heavily confused by 
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hyperfine structure, see iLiszt k Lucas! (|2001h for a sur- 
vey discussion of the CN-bearing species. The features at 
— 17.5 km s and — 10.5kms _1 are more chemically com- 
plex, with much more readily detectable amounts of the 
CN-family and such sulfur-bear ing species as CS, SO , H 2 S 
and HCS+ (ILiszt k Lucas]l2002h or H9CO and NH-, dNash 
19901: iMarscher et alJll99alLiszt k LucaJl995llLiszt et al 



2006). 



The top panel of Fig. Q] shows 1' HPBW CO and 
13 CO J=l-0 emission spectra from the former NRAO KP 
12m telescope, which we used (in conjunction with the 
absorption spectrum) to check the CO column density 
and excitation tempe ratures toward this and other sources 
(jLiszt k Lucasl Il998l ) . The two rather chemically-different 
absorption line components at —10.5 and — S.Skms^ 1 are 
heavily blended at 1' resolution, unlike the clearly double 
structure of the CO and other absorption. This is atypi- 
cal of our comparisons in many other lines of sight between 
the absorption spectra of different species and the CO emis- 
sion profiles when the latter are so strong. It is also curious 
because, being chemically so different, it would seem rea- 
sonable to believe that the two kinematic components were 
unrelated physically along the line of sight. Alternatively 
the —10.5 and — 8.5 km s" 1 components could be related 
and mixed with or enveloped by warmer optically thin gas 
which fills in the emission spectrum at — lOkms -1 but not 
the absorption. However, as will be shown here, the two fea- 
tures are indeed part of the same structure and disparity 
between the emission and absorption profiles arises solely 
as the result of beam-smearing at 1' resolution. 

Several of the features (at —14 and — 4.5/ — S.Skms^ 1 ) 
seen so strongly in CO and HCO + absorption are almost 
absent in CO emission. Can these particular clouds re- 
ally be dense enough to form HCO + and the hydrocar- 
bons C2H and 3h2, along with relatively small amounts of 
HNC and HCN, but not d ense enough to fo rm detectable 
amounts of NH 3 or H 2 CO (jLiszt et al.l l2006) or even dense 
enough to excite CO to detectable levels? The clouds at 
— 14, —4.5/ — 3.5 km s -1 which are so weak in CO emission 
have non-negligible J=l-0 optical depths 0.3 - 0.6 in 12 CO, 
implying that their excitation temperatures are less than 
1 K above the cosmic microwave background. Comparably 
low excitation temperatures have recently been found to 
be common in CO seen in the uv absorption spectra of 
bright stars when N(CO) < 2 x lQ 1 ^ cm" 2 and N(H 2 ) 
< 5 x 10 20 cm -2 (|Burgh et all 120071 ; iSonnentrucker et all 
I2007IL the latter equivalent to N(HCO+) <, 10 12 cm" 2 . 
Although one might also wonder if such H2-bearing, CO- 
silent clouds are just very cold, perhaps part of the dark 
molecular baryonic matter whi ch has been hypothesiz ed to 
exist in the disks of galaxies (Pfcnnig er et al.| [l994). the 
H 2 seen in uv absorption along comparable sightlines host- 
ing weakly-excited CO indicates kinetic temperatures of 40- 
80 K and the kinetic temperature inferred from C 2 is only 
slightly lower (ibid). 

In summary, the line of sight towards NRAO 150 is very 
heavily extincted, having an exceptionally large column of 
intervening neutral gas compared to other sources we have 
studied. Nonetheless, the occulting material is not man- 
ifestly dark in a local sense. About 6 mag of visual ex- 
tinction comes from low density (1cm -3 ) gas along the 4- 
5kpc length of the sightline. Comparable intervening col- 
umn density of H-nuclei comes from 5-6 diffuse clouds 



(A v ~ 1 mag each) , which display a surprisingly complex 
chemistry relative to their modest inferred column density. 

In an attempt to characterize more fully the nature 
of the host gas clouds, we undertook to map CO emis- 
sion around the NRAO150 sightline. However, targeting 
NRAO 150 involves making one serious tradeoff: in return 
for access to kinematic complexity along with the ability to 
map several strongly-absorbing clouds simultaneously and 
the chance to address various puzzles raised by our previous 
observations, we forfeit precise knowledge of the distances 
of the intervening clouds. Only an angular scale may read- 
ily be attached to our present conclusions, a fact typical of 
observations in the galactic plane. For reference, we note 
that the line of sight velocity gradient in the disk, assum- 
ing a flat rotation curve with Rq = 8.5kpc and a rotation 
velocity of 220 km s -1 implies 1) a line of sight distance of 
2.5kpc for a velocity of — 20kms~ 1 (i.e. at the blue edge 
of our strong absorption features) and 2) a distance uncer- 
tainty of ±0.75kpc associated with the typical cloud-cloud 
velocity dispersion of 6kms~ 1 . At a distance of lkpc, 1' 
corresponds to 0.3pc. 

3. Observations and data reduction 

3.1. Observations 

We used the IRAM Plateau de Bure Interferometer (PdBI) 
to achieve higher spatial resolution in mapping CO emission 
from the gas toward NRAO150. However, an interferometer 
filters out low spatial frequencies of the source distribution 
and to supply the missing but much-needed information we 
observed the same region of the sky with the IRAM-30m 
telescope. In this Section we describe in detail how these 
observations were taken, reduced and merged to form the 
on-the-fly (OTF) and hybrid synthesis results. 

3.1.1. IRAM Plateau de Bure Interferometer observations 

PdBI observations dedicated to this project were car- 
ried out with 5 antennae in the CD configuration (base- 
line lengths from 24 m to 229 m) from July to December 
2005. Three correlator bands of 20 MHz were concatenated 
around the 12 CO J=l— frequency to cover a ~ 50kms _1 
bandwidth at a resolution of ~ O.lkms -1 . Two addi- 
tional correlator bands of 320 MHz were used to measure 
the 2.6 mm continuum over the 580 MHz instantaneous IF- 
bandwidth available with this generation of receivers. 

We observed a seven-field mosaic. The mosaic was cen- 
tered on the position of NRAO150 and the field positions 
followed a compact hexagonal pattern to ensure Nyquist 
sampling in all directions. This mosaic was observed for 
about 24 hours of telescope time with 4 to 6 antennas. 
Taking into account the time for calibration and after filter- 
ing out data with rms tracking error larger than 1" (com- 
puted over the scan length), this translated into on-source 
integration time of useful data of 3.3 hours in the D con- 
figuration and 7 hours in the C configuration for a full 6- 
antenna array. Typical 2.6 mm resolution is 5.8" in the D 
configuration and 2.6" in the C configuration. 

3.1.2. IRAM-30m singledish observations 

The 12 CO J =1-0 and J=2-l lines were observed simul- 
taneously during about 4 hours of good summer weather 
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Table 1. Observation parameters. The projection center of all the data matches the position of NRAO150: 0:2000 = 
03 h 59 m 29.747 s , £ 20 oo = 50°57'50.16". 



Molecule 


Transition 


Frequency 


Instrument 


Config. 


Beam 


PA 


Vel. Resol. 


Int. Time 


T S y S 


Noise" 


Obs. date 






GHz 






arcsec 


O 


km s — 1 


hours 


K 


K 






J=l-0 


115.271202 


PdBI 


D 


6.0 X 5.5 


82 


0.2 


3.3/12.5 


230 


0.20 


Jul. 2005 


12 CO 


J=l-U 


115.271202 


PdBI 


C 


2.7 x 2.5 


12 


0.2 


7.0/12.0 


230 


0.52 


Dec. 2005 



a The noise values quoted here are the noises at the mosaic phase center (Mosaic noise is inhomogeneous due to primary beam 
correction; it steeply increases at the mosaic edges). 
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Frequency 
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# Pix. 


FcS 


B cB 
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Resol. 


Int. Time" 


J- ays 


Noise 5 


Obs. date 






GHz 










arcsec 


km s 


hours 


K 


K 




LA CO 


,7=1-0 


115.271202 


30m/AB100 


2 


0.95 


0.75 


22.5 


0.20 


2.4/3.6 


240 


0.26 


Sep. 2005 


12 CO 


J=2-l 


230.538000 


30m/AB230 


2 


0.91 


0.52 


22.5 


0.20 


2.4/3.6 


320 


0.16 


Sep. 2005 



a Two values are given for the integration time: the on-source time and the telescope time. 
b Noise values estimated at the position of NRAO150. 



Table 2. NRAO150 flux in Jy. 



Date 


2.6 mm 


1.3 mm 


17.07.2005 


2.8 


1.6 


28.07.2005 


2.8 


1.6 


01.08.2005 


2.9 


1.7 


12.12.2005 


2.4 


1.4 



(~ 4 mm of water vapor) in September 2005. The four avail- 
able single-side band mixers were used (two per frequencies: 
AB100 and AB230) in the OTF observing mode to map 
a field-of-view of 220" x 220" centered on the position of 
NRAO150. The dump time was 1 second and the scan- 
ning speed ~ 3"/sec. The field-of-view was covered twice 
by raster-scanning along the RA axis and once by scanning 
along the Declination axis. The HPBW of the telescope for 
the J=l-0 line is 22" and the rasters were separated by 10", 
providing Nyquist sampling at 2.6 mm but undcrsampling 
at 1.3 mm. We estimate the 30m position accuracy to be 
~ 2". 

We used the ON-OFF switching mode and separately 
observed the reference position at (+400", 0") during 62 
minutes using frequency switching mode with a frequency 
throw of 7.7 MHz (optimal to suppress standing waves) to 
check for the presence of signal. Emission was detected at 
A2.6mm and later added back into the on-the-fly emission 
profiles. 

3.2. Data reduction 

The data proces sing was done with the GILDAsQ software 
suite (|Petvll2005h . 

3.2.1. Singledish calibration 

The IRAM-30m data were processed inside GILDAS/CLASS 
software. They were first ca librated to the Tt scale using 
the chopper wheel method (jPenzias fc Burrud[l973l ). and 
finally converted to main beam temperatures (T m b) using 
the forward and main beam efficiencies (F c g & i? c ff) dis- 
played in Table [TJ The resulting amplitude accuracy is ~ 

1 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS softwares. 



10%. The line window was fixed between —25 and 10 km s _1 
and a second order baseline was subtracted from each on- 
thc-fly spectrum. The off-position frequency-switched spec- 
trum was fitted by a combination of Gaussians after base- 
line subtraction and this fit was then added to every on- 
the-fly spectrum. 

To compensate for the presence of absorption in the 
OTF-data near the continuum background source, the 
absorption spectra m easured several years ago at PdBI 
(jLiszt fc Lucas!ll998D were scaled to the flux of NRAO150 
at the observation date (Table [2] here), converted to main 
beam temperature, convolved by a Gaussian whose full 
width at half maximum corresponds to the natural res- 
olution of the 30m at the observing frequency, and (at 
last) subtracted from each on-the-fly spectrum. The result- 
ing spectra were finally gridded through convolution by a 
Gaussian. The 12 CO J=2-l data were gridded at the same 
spatial resolution (e.g. 22.5") as the 12 CO J =1-0 data on 
account of the observational undcrsampling at the higher 
frequency 

3.2.2. Inteferometric calibration 

The first steps of the standard calibration methods imple- 
mented inside the GILDAS/CLIC software were used for the 
PdBI data. The radio-frequency bandpass was calibrated 
using observations of 3C354.4 leading to an excellent band- 
pass accuracy because its 2.6 mm flux was larger than 20 
Jy at the epochs of the observations. The flux of NRAO150 
was determined against the primary flux calibrator used at 
PdBI, i.e. MWC349. The resulting flux accuracy is - 15%. 

We then took advantage of the presence of the strong 
NRAO150 point source continuum (Table [2]) to get a very 
accurate phase and amplitude calibration, uv tables for the 
2.6 mm continuum and the 12 CO J=l— were created and 
time dependencies of the phase and amplitude of the con- 
tinuum were self-calibrated per baseline through the use 
of the GILDAS/UV_GAIN task. The temporal phase and am- 
plitude gain derived during this self-calibration step were 
applied to the 12 CO J =1-0 uv-tables through the use of 
the GILDAS/UV_CAL task. Those last two steps were done 
for each individual mosaic field taking into account that 
the flux of NRAO150 is attenuated differently by the inter- 
ferometer primary beam depending on the field position. 
They were also done on each individual observation ses- 



J. Pety, R. Lucas, and H. S. Liszt: Bright, turbulent carbon monoxide 



5 



sion to take into account the temporal variations of the 
NRAO150 flux. Finally, the 12 CO J=l-0 absorption spec- 
trum was added back to the data after rescaling to the flux 
of NRAO150 at the epoch of observation. 

3.2.3. Short-spacing processing, imaging and deconvolution 

Following iGueth et all (jl996l ), the GILDAS/MAPPING soft- 
ware and the singledish map from the IRAM-30m were 
used to recreate the short-spacing visibilities not sampled 
at the Plateau de Burc. These were then merged with the 
intcrferometric observations. Each mosaic field was then im- 
aged and a dirty mosaic was built combining those fields in 
the following op timal way in terms of signal-to-noise ra- 
tio (|Guethll200lD 



Bi{a,S) 



Bi(a,5f 



In this equation, J(a, d) is the brightness distribution in 
the dirty mosaic image, Bi are the response functions of the 
primary antenna beams, Fi are the brightness distributions 
of the individual dirty maps and di are the correspond- 
ing noise values. As may be seen in this equation, the dirty 
intensity distribution is corrected for primary beam attenu- 
ation, which make the noise level spatially inhomogeneous. 
In particular, noise strongly increases near the edges of the 
field of view. To limit this effect, both the primary beams 
used in the above formula and the resulting dirty mosaics 
are truncated. The standard level of truncation is set at 
20% of the maximum in MAPPING. The dirty image was de- 
convolved using the standard Hogbom CLEAN algorithm. 
The resulting data cube was then scaled from Jy/beam to 
T m b temperature scale using the synthesized beam size (see 
Table H]). 

Although we also obtained data in the C configuration 
of PdBI, we present here only the combination of the more 
compact D configuration with the IRAM-30m singledish 
data. Indeed, we were not able to correctly deconvolve the 
combination of C + D configuration and singledish data. 
We guess that this is due to the too low signal-to-noisc 
ratio of the C configuration data. 

3.3. Field of view and finally-achieved resolution 

The resolution of the shown OTF singledish maps is 22.5" 
both for the 12 CO J=l-0 and J=2-l lines. As noted in 
Tabic 1 , the spatial resolution of the hybrid synthesis results 
shown here (using the D-array data) is 5.8" (6.0" x 5.5"). 

The 220" extent and 90" diameter of the OTF and hy- 
brid maps correspond in angle to the HPBW of typically- 
illuminated telescopes (HPBW = 1.2A/D) of diameter D 
= 2.9m and 7.2m in the J=l-0 line and the 5.8" synthesis 
beam corresponds to D — 109m. 

3.4. Relative intensities 

The validity of the discussion depends on the consistency 
of the intensity calibration across several instruments and 
over a wide variety of spatial scales. For the Tr* scale em- 
ployed at the 12m telescope, the antenna temperature is 
standardly scaled up by 1/0.85 representing the fraction of 
the forward response subtended by the Moon. To put the 
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Fig. 2. Evolution of J=T0 12 CO brightness with angular 
resolution. At bottom, the profile from the former NRAO 
12m telescope is compared with hybrid synthesis data con- 
volved to 1' HPBW. At top the hybrid data centered on 
the continuum source are displayed at full (5.8") and 30m 
(22.5") resolution. Note that the hybrid synthesis profile 
much more closely resembles the absorption spectra shown 
in Fig. [1] at bottom. 



data on its standard main-beam scale T m b, IRAM-30m an- 
tenna temperatures are scaled up by 1/0.75 (see Table 1 
here) representing the fraction of the forward response in 
the main beam. 

To make the point that these intensity scales are consis- 
tent, Fig. [1] and [5] compare the profiles toward NRAO150 at 
the same angular resolution from different instruments: 1) 
22.5" spectra from the IRAM-30m and the smoothed syn- 
thesis data and 2) 60" spectra from the Kitt Peak 12m and 
the smoothed synthesis data. The thin line (red) spectrum 
in the upper part of Fig. [2] displays the convolution of the 
hybrid synthesis data to the resolution of the IRAM-30m 
telescope. The agreement with the IRAM-30m spectrum 
shown in the top panel of Fig.[T]demonstrates the reliability 
of the algorithm used to merge the singledish and intcrfero- 
metric data. At bottom in Fig. [2l we also show the result of 
convolving the hybrid synthesis data to the 60" resolution 
of the 12m telescope without any further amplitude rescal- 
ing. However coincindentally, the convolved hybrid synthe- 
sis data reproduce the 12m telescope integrated intensity 
to within 4 percent in either strong feature as well as the 
detailed profiles themselves (although there is a very slight 
velocity shift in the gas around — lOkrns -1 ). Therefore we 
are confident that the very strong lines seen here at high 
resolution, 12-13 K, are not artifacts of instrumental cali- 
bration. 

Finally, the 22.5" resolution emission profile from the 
IRAM-30m in the top panel of Fig. Q] is more than twice 
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as strong as that from the 12m (not just 0.85/0.75) and 
that from the hybrid synthesis at the top of Fig. [2] is 30% 
stronger still. This profile evolution arises from increased 
spatial resolution. 

4. New Observational Results 

4.1. Emission spectra toward NRAO150 

The upper part of the top panel of Fig.[T]compares the 22.5" 
resolution IRAM-30m spectrum from the present work with 
its V KP counterpart toward the continuum source. To the 
extent that the —17.5 and — 8.5kms _1 emission features 
resemble themselves as the linear resolution increases by a 
factor of 3, it might be inferred that the CO emission shows 
little structure on sub-arcminute scales. Conversely, the 
brightness increases by almost a factor 2 at — 10.5 km s -1 , 
indicating stronger sub-structure. 

While the 22.5" resolution IRAM-30m spectrum re- 
solves some but not all of the disparity between the 1' KP 
CO emission profile and the distinct doubling of the ab- 
sorption components at —10.5 and — 8.5kms _1 , the 6" on- 
source hybrid synthesis profile in Fig.[2]much more strongly 
resembles the absorption profile around — lOkrns -1 . So, af- 
ter much work, we are finally able to say that the disparity 
between the 1' KP12m profile and the absorption is due 
to spatial gradients and the coincidental positioning of the 
continuum source with respect to the foreground gas. 

4.2. Summary maps of emission integrated over each 
absorption feature 

To understand the overall behavior, Fig. [3] shows summary 
maps integrated over five velocity ranges corresponding to 
the features seen in absorption. 

The feature at — 17.5kms _1 is present over the full ex- 
tent of the region mapped and/or synthesized but with sub- 
stantial intensity contrast, of order 4:1 even in the middle, 
top panel of Fig. [3J At 6" resolution at right, there is emis- 
sion toward and around the continuum and a somewhat 
higher contrast, from a low of 2.1Kkms~ 1 (just north of 
NRAO150) to 14.5 K km s -1 near the southwest rim of the 
synthesized region. 

The weak emission features at —14 and near — 4kms~ 1 
occupy relatively small fractions at the edges of the area 
mapped at 5.8" resolution in Fig. [3] and approach but do not 
extend across the continuum position. Stronger emission 
exists within the IRAM-30m map outside the hybrid syn- 
thesis field of view but the beams employed here are mostly 
unfilled when pointed on-source: there are no isolated hot 
spots within the synthesized region which were lost to beam 
dilution in the singledish data. We can thus deduce that any 
difficulty in detecting emission around NRAO150 in these 
features is purely an excitation effect, because the optical 
depths measured from the absorption profiles are of order 
0.5 in either case. 

As suggested by the comparison between the 30m and 
12m profiles in Fig. [TJ the sharpest spatial structure near 
the continuum source occurs in the panels centered at 
— 10.5kms _1 in Fig. [3j The ridge of emission which passes 
near NRAO150 appears to have been at most barely re- 
solved by the 30m telescope at — 10.5kms _1 , while a very 
similar structure at — 8.5kms _1 in the OTF maps is more 
spatially extended and better resolved. 




-15 -10 
V LSR (KM S" 1 ) 

Fig. 4. Hybrid synthesis CO J=T0 (5.8" resolution) pro- 
files at selected locations of particularly strong emission. 
Top: Profiles selected from a map of peak emission near 
— lOkms^ 1 . Bottom: Selected from a map of peak emis- 
sion near — 17kms _1 . Profiles are labelled by their dis- 
placement (Aa, A5) in arcseconds relative to the contin- 
uum source. 



Perhaps the clearest result in Fig. [3] is the transverse 
resolution and isolation of the ridge at — 10.5kms _1 near 
NRAO150 in the hybrid synthesis data at right, and, to 
a slightly lesser extent, at — 8.5kms _1 as well. Over the 
region of the synthesis the integrated intensity contrast is 
high, above 16:1, but the emission does not break up into 
a patchwork of unresolved and/or isolated bright and dark 
regions. As discussed below, uniformity is to some extent an 
artifact of the high optical depth in the middle of the ridge, 
enhanced by integration over velocity which suppresses spa- 
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Fig. 3. Summary of mapping information of the diffuse gas around NRAO150. Each panel shows the 12 CO emission 
integrated over the velocity range of one absorption feature visible in Fig. [IJ i.e. from top to bottom the —17.5, —14, 
— 10.5, — 8.5kms _1 and, combined, the —4.5 and — 3.5kms _1 features. The position of NRAO150 is marked by a cross 
at the map center. The left and middle columns show respectively the 12 CO J=2-l and J =1-0 emission observed with 
the 30m singledish antenna while the right column shows the 12 CO J=l-0 emission at higher resolution from hybrid 
synthesis data (i.e. the combination of the 30m and PdBI instruments). The 30m 12 CO J =2-1 map has been smoothed 
to 22.5". The field of view of the hybrid synthesis maps is about half that of the singledish data: As a guide, the red circle 
has the same 45"-radius in all panels. The resolution of the hybrid synthesis data (5.8") is almost 4 times better than 
that of the singledish data, 22.5". Contour levels are spaced by 2 Kkms -1 , but for the first contour for the singledish 
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tial structure even at the northern edge of the ridge where 
(see Sect 4.5 and Fig. [7]) it is heavily kinematically struc- 
tured. 

Taken together, similarity in the third and fourth row 
of Fig. [3] strongly suggests that the features seen at —10.5 
and —8.5 kms -1 are actually part of the same structure, 
no matter that they are distinct in absorption, separated 
in emission at sufficiently high spatial resolution and sub- 
stantially different in the richness of their chemistry. 

4.3. Channel maps at 1 km s" 1 resolution 

Fig. [10] and [TT| (available on-line) show a series of channel 
maps with a resolution of lkms -1 for the IRAM-30m and 
PdBI data. The weakly-emitting gas at —14, —4.5 and/or 
—3.5 kms" 1 is only marginally seen near the field center 
and an East- West velocity gradient across NRAO150 is re- 
solved in the —17.5 km s" 1 gas with more positive velocities 
seen to the South and East. The blue wing of emission at 

— 11 kms -1 is disposed on either side of the source (some- 
what filling in the gaps of the distribution at —10 kms" 1 ). 

Increasing the velocity discrimination to this extent still 
does not show substantially smaller scale spatial struc- 
ture in the emission pattern with the possible exception 
of the panel at —Skins" 1 , in the near wing of the strong 
emission feature. The high optical depths measured in the 

— 17.5 kms" 1 feature and in the components at —10.5 and 
—8.5 kms" 1 may impede viewing such structure: No mat- 
ter how sharp the actual internal structure, it can only be 
seen if the line of sight penetrates deeply enough to probe 
it. This behavior is illustrated in greatest detail in Fig. [7] 
and discussed in Sect. 4.5. 



-10.5 & -8.5 km/s -17.5 km/s 
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Fig. 5. Top: Probability distribution function of the line 
brightness in the hybrid synthesis, normalized to unit area. 
Bottom: Fractional area (cumulative distribution) of pix- 
els whose line brightness is higher than the value given on 
the abscissa. Right: for the component at —17.5 kms" 1 . 
Left: for the —10.5 and —8.5 kms" 1 components. The 
histograms have been computed in the line wings (blue 
curves), the line cores (red and green curves) and baseline 
devoid of signal (black curve). For reference, the averaged 
spectra are shown as inset with the same color scheme. The 
vertical dotted lines display the line brightness medians. 



4.4. Bright spots, area filling factors and spatial statistics 

Fig. H] displays profiles selected from hybrid maps of peak 
brightness in the gas around —10 km s" 1 (see the top panel) 
and at —17 kms" 1 . There are some very strong lines in the 
region of the synthesis, T m b= 12-13 K, which was hardly ex- 
pected on the basis of the 3 K lines in the emission spectra 
at 1' in Fig. [TJ As indicated by the spatial offsets labelling 
the spectra in the top panel of Fig. [4) such unusually strong 
lines exist very close to NRAO150, i.e. just outside the on- 
source synthesized beam but still inside the on-source 30m 
beam. These lines imply CO J=l-0 excitation temperatures 
of at least 16-17 K, far stronger than is typical of very dark 
clouds wh ere Tk < 10 K. P erhaps an analogous surprise 
was that of lHeithausenl (|2004f ). whereby a very weak (0.2 K) 
12 CO J=T0 line inadvertently seen at the 30m was subse- 
quently found in high-resolution PdBI maps to consist of a 
much brighter, unresolved spicule. 

The beam efficiency needed to convert observed antenna 
temperatures to actual sky brightness has two aspects; the 
extent of the sky distribution - which antenna lobes are 
subtended - and the degree of uniformity or area filling 
factor within that extent. In fact there is a high degree of 
heterogeneity on all the scales which are resolved in the 
present dataset. Fig. [5] expresses this fact quantitatively by 
showing the fractional area occupied by emission at all ob- 
served levels of brightness. We defined a circle of radius 
45" within the hybrid synthesis map and constructed his- 
tograms of the brightness distribution for all included pixels 
over various channel ranges. At left in Fig. [5] are the differ- 
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Fig. 6. Mean profiles (< T m b >, red lines) and rms spatial 
brightness fluctuation in excess of the random noise rms 



(yRMS(T mb ) 2 - RMS (noise) 2 , shaded histogram). Left: 
For the 5.8" resolution hybrid synthesis data (random noise: 
0.47 K). Right: For the 30m OTF mapping data at 22.5" 
resolution (random noise: 0.33 K). 



cntial and cumulative distributions for the near wings and 
cores of the —10.5 and —8.5 kms -1 features and at right 
are those for the —17.5 kms" 1 feature. 

Cartoons inset in each panel show the profile inter- 
vals which were employed, superposed on the mean profile 
formed over the region. These mean spectra have bright- 
nesses of 3.6 - 4.2 K in the line cores and approximately 
1 K in the selected regions of the line wings. Noise statistics 
are not necessarily the same for channels with and without 
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emission (after deconvolution) but for reference we con- 
structed similar histograms for channels in emission-free 
baseline regions of the spectrum (as shown). 90% of the 
noise occurs below 0.7 K and 95% below IK; the nominal 
single-channel rms noise averaged over the region of the 
synthesis is 0.47 K in emission-free baseline regions of the 
spectrum. 

Pixel histograms for the line cores are strongly popu- 
lated up to 7-8 K and have tails extending to 10-13 K. Even 
pixels in the 1 K line wings have higher-intensity tails ex- 
tending up to 6-7 K, reflecting substantial local excursions 
of the line centroids (see Fig. 2] and [7]). Such wide distribu- 
tion of intensity ensures that fluctuations are a significant 
fraction of the mean. A plot of the true spatial rms varia- 
tion in each channel, i.e. the rms variation above and be- 
yond that attributable to radiometer noise (e.g. 0.47 K av- 
eraged over the field of the the hybrid synthesis and 0.33 K 
in the 30m OTF data), calculated as the square root of the 
quadrature difference between the total and noise rms, is 
shown in Fig. [6] 

For the larger 220" x 220" region of the 30m OTF data 
this excess rms is closely equal to the mean across the entire 
spectrum, except between —9 and — lOkms -1 where the 
mean exceeds the rms by 20-30%. For this gas, increasing 
the spatial resolution to 5.8" yields only a small increase in 
the ratio of mean and excess rms brightness. By contrast, 
the character of the comparison changes drastically for the 
— 17.5 feature at higher resolution, where the excess rms 
changes from 100% to only 50% of the mean at the line core 
when going from the IRAM-30m to the synthesis data. 

The actual pixel statistics of the brightness distribu- 
tions are given in Fig. [5] but fluctuations equal to the mean 
are most commonly recognized as characteristic of Poisson 
statistics, which would be observed from a macroscopically 
transparent ensemble of statistically independent clumps. 
Conversely, large ratios of mean to rms brightness might 
imply that the innate structure has been completely spa- 
tially resolved. However the situation is both more com- 
plex and more interesting, because the synthesis data also 
resolves the kinematic substructure of the features in ways 
which are not reflected in this discussion. 

4.5. A detailed view of line kinematics 

The full extent of the variation in CO emission can only 
be revealed by examining the data at the highest available 
spectral resolution. Fig. [7] explores the kinematics on a se- 
ries of tracks along and across the rather linear structure 
seen around — lOkms -1 . The tracks, shown overlaid on an 
integrated intensity map at top left in the Figure, are sep- 
arated vertically (in Declination) by 4 pixels or 6" which 
is very nearly one synthesized HPBW. Spectra along the 
transverse track p are shown at bottom left, where the ver- 
tical scale is elongated compared to the others at right in 
the figure. At top right is a larger-scale diagram along track 
d using 30m data and below it are shown the hybrid syn- 
thesis data along tracks a-f. Comparing the track at the top 
right in Fig. [7] with those below is a striking deconstruction 
of the singledish lower-resolution data into its interferomet- 
ric high-resolution components. 

The simple, essentially featureless East- West velocity 
gradient in the — 17.5 km s -1 gas at top right in Fig. [7] 
is manifest along strips a-e but the manner in which the 
velocity change occurs is striking at high resolution, with 



abrupt reversals and line-splittings; it is zig-zag and me- 
andering rather than smooth. The same is true for the 
emission around -lOkms , in those strips (a-e) mostly 
above the mid-line of the long trunk. Along track f, the line 
profile over the complex is filled in by stronger and pre- 
sumably very opaque emission seen in the profile labelled 
(-7.5",-13.5") in Fig. H We hypothesize that the line of 
sight penetrates deeply into the medium around —10 km s _1 
at the edge of the trunk along tracks a-d, and that kine- 
matic details largely disappear along the tracks e-f as the 
medium becomes fully opaque in a macroscopic sense; note 
the excursions in the line wing at v < — llkms -1 even 
along tracks e-f. The combination of beam-smearing and 
the strength of emission just south of NRAO150 combine in 
on-source singledish data to obscure the double-lined struc- 
ture apparent in absorption or at high spatial resolution. 

4.6. Rotational excitation 

Insight into the nature of the host gas is hindered by the 
unavailability of high-resolution observations other than of 
12 CO and the lack of high spatial resolution even in the 
J=2-l line. However, in diffuse gas where the CO/H2 ratio 
is small and sharply variable, where the 12 CO/ 13 CO ra- 
tio can be strongly altered by fraction ation, and where the 
rotational excitation is sub-thermal dLiszt fc Lucai 1998; 
I Sonnentrucker et al.ll2007l : iBurgh et al.ll2007HSheffer et al l 
l2007t lLisztll2007al ). many of the underpinnings of typical 
analysis of CO profiles must be questioned in any case. 

Abundance effects can be factored out somewhat by 
studying 12 CO alone, as shown in Fig.[8j Such plots compar- 
ing the J=2-l and J=l-0 intensities confirm what is already 
evident in Fig. [3J that the two distributions are remark- 
ably similar, especially given the sub-thermal excitation in 
diffuse gas. The data indicate a nearly universal ratio of 
intensities, approximately 0.725:1 with relatively small dif- 
ferences among the various features. The — 14kms _1 fea- 
ture, whose excitation is too weak to produce detectable 
emission toward NRAO150 has the same slope in the dia- 
gram as much more strongly emitting features, and only the 
gas at — 4kms _1 seems to have a notably different, smaller 
slope. A near-constant line brightness ratio, with a slightly 
smaller slope, was noted prev iously in low column density 
gas bv iFalgarone et aD (|l998h . see their Fig. 11 and refer- 
ences. The easiest way to explain this ratio is to vary the 
degree of beam-filling of optically thick LTE emission at 
Tk = 14 K. We however explain in the next section why 
this scenario is unlikely here. 



5. Modelling 

5.1. Density 

To get some feeling for the observed properties of the 
CO emission and the density in the host gas, we per- 
formed non-LTE microturbulcnt radiative transfer calcu- 
lations in stationnary models of diffuse cloud. These mod- 
els self-consistently determine the kinetic temperature, H2 
and CO abundances and CO excitation temperature in 
uniform density (n(H) — n(H I) + 2n(H2)) gas spheres. 
The column density thus varies across the sphere, imply- 
ing that the modelled profiles samples very different col- 
umn densities, excitation conditions and c hemical behav- 
iors. These models are described in details in Liszt & Lucas 
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Fig. 7. Position-velocity diagrams across and along the main trunk feature at a position angle of —50°. Top, left: 
Integrated intensity map of 5.8"-resolution hybrid synthesis data summed over the disjoint velocity ranges [—19.. — 
15]kms _1 and [—12.5.. — 7.5] kms" 1 . Contours are shown at levels 5,10,15,..35Kkms -1 . The position of NRAO150 is 
marked by a cross at the map center. Bottom, left: A position-velocity diagram along the perpendicular strip p 
through the position of NRAO150 at 5.8" resolution. Contours are shown at levels 0.5,1,2,4,6,8,10,12 K. Top, right: 
Longer position-velocity diagram at 22.5" resolution straddling strip d, made from 30m data. Contours are shown at levels 
0.5,1,2,4,6,8,10,12 Kkms" 1 . Bottom, right: Position-velocity maps at 5.8" resolution along strips a-f. Coordinates are 
displacement in right ascension relative to NRAO150. Strip d precisely crosses the continuum source. 



(|2000l ) and[Liszt| (|2007bh . They were used bv lLisztl (|2007aD 
to interpret optical observations of CO excitation in dif- 
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Results for the models are given in Fig. [9] and they 
show that the observed line ratios are as expected for non- 
LTE CO emission from cool media of low-moderate num- 
ber and column density n(H) — 64 — 512 cm" 3 , iV(H) < 
3.8 x 10 21 cm" 2 or Ay < 2 mag. The pressures in these 
models are p/k — (2 — 8) x 10 3 K cm" 3 and the H2/H ratio 
varies considerably (increasing with n(H) or with N(CO), 
hence Tb and the other results are line-of sight averages 
like the data). The observational noise levels are sufficient 



to obscure the small differences in line ratio which would 
distinguish between model densities for J=l-0 brightnesses 
below 2 K and useful differences between the model results 
at various total densities are really only apparent in some- 
what stronger lines. 

The gas around — 4kms~ 1 is distinguished in the right- 
most panel of Fig. [5] by its locus below the indicated line of 
slope 0.725 and this deviation is readily interpreted in terms 
of the model results for somewhat lower density perhaps 
n(H) = 64 — 128 cm -3 , although small line ratios also arise 
in the outer regions of model spheres with higher density if 
the column density is small enough. By contrast, the loci for 
the gas at — 17.5 kms -1 and — 10.5 kms are slightly above 
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-19 < v < -16 km/s -15 < v < -13 km/s -13 < v < -9.5 km/s -9.5 < v < -6.5 km/s -6 < v < -3 km/s 
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Fig. 8. Joint distributions of the emission of 12 CO J=2— 1 and J =1-0. The value at a given position of these joint 
histograms is the number of pixels of the input images whose intensities lye in the respective vertical and horizontal bins. 
Contour levels are set to 1, 2, 4, 8, 16, ... 512 points per pixel. The joint distributions have been computed and plotted for 
each velocity components. A (red) line of slope of 0.725 is shown for reference. Only the emission around v = — 4kms~ 1 
in the rightmost panel departs noticeably from the mean behavior. 
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Fig. 9. Non-LTE models of CO brightness. The kinetic tem- 
perature, molecular hydrogen and carbon monoxide abun- 
dance are self-consistently determined for gas spheres of 
constant total density n(H) = n(H I) + 2 n(H2) = 64 .. 
512 cm -3 as labelled. The CO line brightnesses are found 
from microturbulent radiative transfer across the faces of 
individual models. Each series of points shows the result 
for a different central column density varying in steps of 2 
over the range 2 x 10 20 cm~ 2 < 2V(H) < 4 x 10 21 cm~ 2 . A 
slope of 0.725 is shown by a dashed (red) line. 



the line at intermediate brightness, arched like the model 
results, suggesting density as high as n(H) i; 256 cm~ 3 . 
Even these are very low densities by the standards of molec- 
ular emission-line work, but they are at the high end of the 
range for diffuse clouds with a substantial residual fraction 
of atomic hydrogen. Somewhat lower d ensities would have 
been inferred if the collision rates of [Balakrish nan et all 
(2002) for CO excitation by H-atoms were employed, but 
these seem now to have been discredited (IShepler et al.l 
I2007IL 



5.2. Bright peaks in diffuse gas 



In emission, any single-dish profile can only be regarded 
as a trade off between the brightness and the beam fill- 
ing factor, only the product being actually observed. On 
the other hand, in absorption, once the apparent opti- 
cal depth is large, the filling factor must approach unity. 
Because the —17.5 and — 10.5 km s -1 features have large 
apparent optical depths in absorption toward NRAO150 
where they are somewhat weaker in emission, they must 
occult the continuum fully (at —17.5 kms" 1 ) or very nearly 
so (— 10.5kms _1 ). In this case, it is likely that the synthe- 
sized beam area is filled, because, if the medium were very 
patchy there would be a high chance that a background 
source could escape occultation and absorption even when 
a foreground cloud was visible in CO emission. Because this 
does not occur in our work we infer that the beams used 
to probe the CO are filled or very nearly filled by the op- 
tically thick emission, leaving little room for an increase in 
the beam filling factor. 

Assuming that the beam is filled, the models illustrated 
in Fig. [S] shows that 10 K and brighter CO J=l-0 lines 
may be produced in diffuse clouds. These occur for n(H) <; 
256 cm- 3 , A(H 2 ) £ 10 21 cm~ 2 N{CO) £ 10 16 cm- 2 when 
the free gas phase carbon is on the verge of recombin- 
ing from C + to CO and the fractional abundance X(CO) 
varies very rapidly with either the number or column den - 
sity of hydrogen, see Fig. 1 and Fig. 6 of iLisztl (|2007aD . 
Even though the CO lines are already optically thick, ob- 
servations of CO in absorption in diffuse gas show that the 
integrated brightness of the J = 1-0 line is proportional 
to the CO column density (ibid). This observational fact 
is inde ed typical of subthermal e xcitation as predicted long 
ago bv lGoldreich fc Kwanl (| 19741 ). Hence, a small increase in 
number or column density of hydrogen at the above thresh- 
old may greatly increase the CO column density, which in 
turn produces brighter optically thick CO lines. 

Bright CO lines around NRAO150 thus most probably 
reflect the variation of a chemical process, i.e. the C + -CO 
conversion. However, the ultimate cause of the variations of 
this chemical process in such a limited field of view remains 
uncertain. 
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6. Summary; prospects for future work 

After a long series of investigations into the absorption 
line chemistry of galactic diffuse gas coincidentally seen 
toward distant radio point sources at low-moderate galac- 
tic latitude, we undertook to image the absorbing medium 
by mapping 12 CO in emission around one source at the 
highest practicably obtainable spatial resolution. There are 
many questions which might be asked of such data and the 
comparisons between relatively microscopic (sub-arcsecond 
pencil absorption beams and the much larger (6 — 60") fan 
beams in emission. The present work is a small step in this 
direction. 

Some of the results are rather direct, for instance con- 
cerning the various disparities and differences between the 
absorption profiles and the existing 60"-resolution KP CO 
data (Fig. [T|). An atypical disparity between the strong 
emission and absorption profiles at — lOkms -1 was shown 
to arise because of insufficient spatial resolution in the 
emission and not, for instance, due to the presence of a 
warmer (more highly-excited) CO component seen in emis- 
sion with low optical depth. The weakness of the —14 and 
—4.5/ — 3.5 kins" 1 emission lines toward NRAO150, which 
are seen in absorption at moderate optical depth (0.5) in 
Fig. [U arises because the regions of sufficient collisional ex- 
citation and stronger emission are somewhat removed from 
the continuum position. These clouds, which are common 
in our previous absorption surveys, are not a previously- 
hidden large scale molecular component revealed in absorp- 
tion, but they do highlight the robustness of the chemistry 
into regions where the density and column density are too 
low to produce much rotational excitation, even in CO. This 
is in fact the norm in optical absorption studies of CO. 

We explored the spatial structure of the gas in several 
ways, i.e. by examining the maps at 22.5" and 6" resolu- 
tion (Fig. © and by calculating the spatial statistics both 
as profile rms over the mean velocity profile (Fig. [6]) and as 
channel by channel probability distributions over the syn- 
thesized region (Fig. [5]). At the lower resolution, the ve- 
locity profile of rms fluctuations due to spatial structure, 
calculated as the channel by channel quadrature difference 
between the total rms and that seen in line-free channels, 
closely resembles the mean profile, except in a narrow in- 
terval near the core of the — lOkms -1 line. Such statis- 
tics are quasi-Poisson in character, consistent with a ran- 
dom macroscopically transparent ensemble of independent 
clumps. However, when the gas is studied at higher res- 
olution, where there is an even wider range of observed 
brightness at all velocities, the statistics in the core of the 
— 17.5kms _1 noticeably change character and the rms is 
only 50% of the mean (Fig. |BJ>. 

In this work, it appears that every profile seen at 22.5" 
with a linewidth greater than that attributable to the local 
sound speed (full-width at half-maximum = 0.96 km s -1 
for a pure-H2 gas with Tk = 40 K) is found to be com- 
posite when viewed at higher resolution: decomposable 
into multiple narrower components and/or having spatially- 
resolved velocity gradients in narrower-lined gas. Whether 
this would continue in the broadest components seen at 
5.8" resolution is an open question. 

The full extent of the emission structure is realized only 
when the highest spatial and velocity resolution are viewed 
together in position-velocity tracks along portions of the 
gas where the medium may be supposed to be macroscop- 



ically optically thin (Fig. [7]) . A series of tracks parallel to 
but somewhat displaced from the peak of a particularly well 
defined filamentary or trunk structure in the gas around 

— 10 kms -1 was used to show the substructure in what are 
almost featureless position- velocity diagrams at 22.5" res- 
olution. The high resolution version resolves a monotonic 
and rather bland velocity gradient in the —17.5 kms -1 fea- 
ture into a striking series of kinks and abrupt reversals 
which strongly re semble simulations of turbulent flow of 
the sort shown by iFalgarone etai] (fl99l . 

The flat profile core of the —10 kms -1 component seen 
at 60" resolution, which is slightly split at 22.5" resolu- 
tion (Fig. [1] and Fig. [5]), shows an underlying triple nature 
(Fig. 2] and Fig. [7]) whose individual components also show 
this striking pattern of velocity reversals when the medium 
is macroscopically thin and show some excursions in the 
line wing where it is thick. In this connection, note that, 
in our data, it was not necessary to avoid the strongly- 
emitting optically thick line cores (at least they are opti- 
cally thick toward NRAO150) in order to see such struc- 
ture. Furthermore, as noted just before, when studied at 
lower angular resolution, regions of the line core exhibit 
the same fractional level of fluctuation as the line wings for 
the component at —17.5 kms -1 . At — lOkms -1 the level of 
fluctuation is not very different at the two resolutions, and 
is still relatively large in the line core. 

In the large, the CO excitation in the host gas sampled 
in the J=l-0 and J=2-l lines at 22.5" resolution (Fig. [5]) is 
readily explained as arising in low-moderate density n(H) 
= 50-250 cm" 3 diffuse gas (Fig.©. The 12-13 K 12 CO lines 
seen in Fig. |H at the high end of the brightness distribution 
(Fig.© are probably regions of higher density (n(H) = 300 - 
500 cm" 3 ) and column density N(H 2 ) > I0 21 cm" 2 , N(CO) 
> 10 16 cm~ 2 where C + -CO conversion is more advanced 
even while the temperature is still elevated, Tk >> 10 K. 
It is something of a coincidence that the 12 CO brightnesses 
observed in diffuse and dark gas are so similar, given the 
profound differences in abundance and excitation: excita- 
tion in diffuse clouds is considerably sub-thermal and the 
temperatures are high, typically above 30 K, and only a 
small fraction of the gas-phase carbon nuclei reside in CO 

- typically 5% or less. This coincidence will complicate the 
interpretation of CO emission in general. 

While absorption studies enable very precise determi- 
nation of various chemical tracer column densitites, they 
can currently be done only on a limited number of pen- 
cil beam lines of sight. Those studies thus miss impor- 
tant information about the environment such as the geom- 
etry and the kinematics of the host gas, information which 
could help to unravel remaining m ysteries like the abun- 

e.q. Falgaronc c t all 120061: 



dance of HCO+ or CH+ (see 



Iffilv-Blant fc Falgarondl2007t iHilv-Blant et al.ll2008ft . This 
paper is the first of a series in which we image diffuse gas 
around lines of sight which have previously been studied in 
absorption. In this series our observations will range from 
the highest possible angular resolution in small fields of 
view (as as present) to lower angular resolution in fields of 
view wide enough to contain the whole structure of the host 
gas. In this way we will relate such different gas and dust 
tracers as H I, CO, HCO + , extinction, etc. using whatever 
means are available to render the diffuse absorbing material 
both more readily visible and more completely understood. 
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Fig. 10. Finer-scale singledish maps of integrated 12 CO brightness. The 30m on-the-fly mapping observations used as 
short-spacing input to the hybrid synthesis are summed over the 5-channel (1.015 kms -1 ) intervals most closely centered 
about the velocity indicated at the upper left in each panel and the integrals are scaled down by 1/1.015 to reflect mean 
brightness over a 1 kms -wide interval. Contours are shown at levels 1,2,4,6,8, lOKkms -1 . The position of NRAO150 
is marked by a cross at the map center and the region synthesized is shown as a circle. Note the velocity gradient in the 
— 17kms _1 feature across the position of NRAO150. The 22.5" HPBW is shown inset at the lower left in each panel. 
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Fig. 11. Finer-scale hybrid synthesis maps of integrated 12 CO brightness. The 5.8" hybrid synthesis profiles are summed 
over the 5-channel (1.05 kms -1 ) intervals most closely centered about the velocity indicated at the upper left in each 
panel and the integrals are scaled by 1/1.05 to reflect mean brightness over a 1 kms -1 -wide interval. Contours are shown 
at levels 1,2,4,6 and 8Kkms _1 . The position of NRAO150 is marked by a cross at the map center. The 5.8" synthesized 
HPBW is shown inset at the lower left in each panel. 
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Fig. 12. Left: Same as Fig. [5] but using the 30m data only. Right: Same as right but using the total field-of-view of the 
30m data, i.e. 200" x 200". 
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Fig. 13. Moments of the components — 17.5 km s" 1 and —10.5 and — 8.5kms _1 together observed at two different res- 
olution. From top to bottom: Peak temperature, Integrated area, Centroid velocity, 2nd order moment (proportional 
to the FWHM for a Gaussian line) and the modulus of the gradient of centroid velocity (for a flat rotation curve of 
the Galaxy, — 20kms~ 1 corresponds to 2.5kpc in the direction of NRAO150 and at 2.5kpc, 0.1 km s^ 1 /" corresponds to 
8 kins" 1 pc _1 ). 
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Fig. 14. Maps of spectra around the profiles shown in Fig. Q] sorted by distance to NRAO150 from top to bottom. 
Left: The -17.5 kms" 1 feature shown around (+3", -3"), (+10.5", -10.5") and (-7.5", -13.5"). Right: The -10.5 
and -8. Skins" 1 features shown around (-27", -27"), -28.5", +25.5" and (+36", -24"). 



